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Abstract 

Background: The mismatch repair (MMR) pathway plays an important role in the maintenance of the genome 
integrity, meiotic recombination and gametogenesis. This study investigated whether genetic variations in MMR 
genes are associated with an increased risk of sperm DNA damage and male infertility. 

Methods: We selected and genotyped 21 tagging single nucleotide polymorphisms (SNPs) in five MMR genes 
{MLH1, MLH3, PMS2, MSH4 and MSH5) using the SNPstream 12-plex platform in a case-control study of 1,292 
idiopathic infertility patients and 480 fertile controls in a Chinese population. Sperm DNA damage levels were 
detected with the Tdt-mediated dUTP nick end labelling (TUNEL) assay in 450 cases. Fluorescence resonance 
energy transfer (FRET) and co-immunoprecipitation techniques were employed to determine the effects of 
functional variants. 

Results: One intronic SNP in MLH1 (rs4647269) and two non-synonymous SNPs in PMS2 (rsl 059060, Ser775Asn) 
and MSH5 (rs2075789, Pro29Ser) seem to be risk factors for the development of azoospermia or oligozoospermia. 
Meanwhile, we also identified a possible contribution of PMS2 rsl 059060 to the risk of male infertility with normal 
sperm count. Among patients with normal sperm count, MLH1 rs4647269 and PMS2 rsl 059060 were associated 
with increased sperm DNA damage. Functional analysis revealed that the PMS2 rsl 059060 can affect the 
interactions between MLH1 and PMS2. 

Conclusions: Our results provide evidence supporting the involvement of genetic polymorphisms in MMR genes 
in the aetiology of male infertility. 



Background 

Infertility remains a major clinical problem that occurs 
in 10 to 15% of couples worldwide [1], and male factor 
infertility accounts for 40 to 50% of all infertility cases 
[2]. Although several causes have been identified for 
impaired male fertility [3], the aetiology remains 
unknown in nearly half of all cases. Currently, a large 
amount of attention is being paid to the potential effects 
of sperm DNA damage on male infertility [4]. DNA 
damage in the male germ line appears as a risk factor 
for adverse clinical outcomes, including poor semen 
quality, low fertilization rates, impaired pre-implantation 
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development, miscarriage and an increased risk of mor- 
bidity in the offspring [5-7]. 

Although the clinical significance of testing sperm 
DNA integrity has been clearly emphasized, the origin 
of DNA damage in spermatozoa is poorly understood. 
One mechanism is that deficits in the DNA repair sys- 
tem during spermatogenesis can have negative effects 
on the integrity of sperm DNA [8,9]. Our previous data 
have provided strong evidence that some genetic poly- 
morphisms in genes involved in DNA repair were asso- 
ciated with the development of sperm DNA damage and 
male infertility [10-13]. 

Among all DNA repair mechanisms, DNA mismatch 
repair (MMR) plays a critical role in the maintenance of 
genetic integrity and malfunctions can lead to various 
cancers in mammals [14-16]. Studies of gene knockout 
mice indicate that several members of the MMR family 
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also participate in the meiotic recombination process 
and are involved in gametogenesis [17,18]. Three MutL 
homologues (MLH1, MLH3 and PMS2) and two MutS 
homologues (MSH4 and MSH5) are involved in this 
process. 

Based on their important physiological functions, these 
five MMR genes are good candidate genes for explaining 
male infertility. Recently, analysis of polymorphic mar- 
kers in candidate genes helped us to understand the 
etiology and the susceptibility of male infertility [19-21]. 
The purpose of this work is three-fold: (1) to examine 
whether MMR gene polymorphisms are associated with 
increased risk of azoospermia or oligozoospermia, (2) to 
ascertain whether genetic variants in MMR genes result 
in sperm DNA damage and, thereby, increase male 
infertility, and (3) to investigate the biological activity of 
the significant functional variants. 

Methods 

Subjects and sample collection 

The study was approved by the Ethics Review Board of 
the Nanjing Medical University. All the studies involving 
human subjects were conducted in full compliance with 
government policies and the Declaration of Helsinki. A 
total of 1,657 infertile patients, diagnosed with unex- 
plained male factor infertility, were drawn from the 
Centre of Clinical Reproductive Medicine between April 
2005 and March 2009 (NJMU Infertile Study). All parti- 
cipants completed an informed consent and a question- 
naire, including detailed information, such as age, 
cigarette smoking, alcohol drinking, tea and vitamin 
consumption, and abstinence time. All patients under- 
went at least two semen analyses, and those with a his- 
tory of orchitis, obstruction of the vas deferens, 
chromosomal abnormalities, or micro-deletions of the 
azoospermia factor region on the Y chromosome were 
excluded [22]. In the final analysis, 1,292 idiopathic 
infertility patients aged 24 to 42 years old were included, 
and were divided into three subgroups: 268 infertility 
patients with non-obstructive azoospermia, 256 inferti- 
lity patients with oligozoospermia (sperm counts < 20 x 
10 /ml) and 768 infertility patients with normal count 
(sperm counts > 20 x 10 /ml). 

The control group included 480 fertile men ranging 
from 25 to 40 years of age who had fathered at least 
one child without assisted reproductive technologies and 
had normal semen parameters. The semen analysis for 
sperm concentration, motility and morphology was per- 
formed following the World Health Organization cri- 
teria [23]. 

SNP selection and genotyping 

We selected the tagging SNPs by using genotype data 
obtained from unrelated Han Chinese individuals from 



Beijing in the HapMap project (HapMap Data Rel 24/ 
Phase II Nov08, on NCBI B36 assembly, dbSNP bl26). 
To examine the gene extensively, we searched the MMR 
genes, including 2,000 bp of the flanking regions both 
upstream and downstream of the gene, using the pair- 
wise option of the Haploview 4.0 software (Mark Daly's 
Lab, Broad Institute of MIT and Harvard, Cambridge, 
Massachusetts, USA). The tagging SNPs were selected 
on the basis of pairwise linkage disequilibrium with a r 
threshold of 0.8 and minor allele frequency > 0.05 to 
capture all the common SNPs. In total, 19 SNPs were 
chosen in these 5 genes. In addition, a non-synonymous 
SNP (rsl799977) in MLH1 and a non-synonymous SNP 
(rs2075789) in MSH5 that cause missense mutations 
were included. 

Genotyping was performed using TaqMan 7900HT 
Sequence Detection System and GenomeLab SNPstream 
high-throughput 12-plex genotyping platform (Beckman 
Coulter, Fullerton, CA, USA). Sequences of forward, 
reverse and extension primers are listed in Additional 
file 1 (Table SI). For quality control, the genotyping was 
done without knowledge of case/control status of the 
subjects, and a random 5% of cases and controls were 
genotyped twice by different individuals, and the repro- 
ducibility was 100%. To confirm the genotyping results, 
selected PCR-amplified DNA samples (n = 2, for each 
genotype) were examined by DNA sequencing and the 
results were also consistent. 

DNA fragmentation analysis 

After a period of 48 to 72 h of sexual abstinence, semen 
samples were collected by masturbation into wide- 
mouthed sterile containers and were delivered to the 
laboratory within 1 h of ejaculation. The diluted samples 
were cooled gradually at 5°C for 2 h, frozen at -70°C for 
Tdt-mediated dUTP nick end labelling (TUNEL) evalua- 
tion. A detailed protocol of the TUNEL assay for 
human sperm has been described previously [24]. 
TUNEL labeling was carried out using a Cell Death 
Detection kit (APO-DIRECT kit; BD Biosciences Phar- 
Mingen, San Diego, CA, USA) according to the manu- 
facturer's instructions. Briefly, semen samples were 
thawed in a 37°C water bath and immediately diluted 
with buffer (0.15 M NaCl, 0.01 M Tris, 0.001 M EDTA, 
pH 7.4) to obtain a sperm concentration of 1 to 2 x 
107ml. Washed sperm was resuspended in 2% parafor- 
maldehyde for 30 minutes at room temperature. After 
rinsing in PBS, samples were resuspended in permeabili- 
zation solution (0.2% Triton X-100, 0.1% sodium citrate) 
for 10 minutes on wet ice. TUNEL reagent (50 ul) was 
added to each sample. For each batch, a negative control 
lacking the terminal deoxynucleotidyl transferase and a 
positive control treated with DNase I were included to 
ensure assay specificity. After incubation for 1 h at 37°C, 
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samples were analyzed immediately by flow cytometry 
(FACSCalibur; BD Biosciences Pharmingen, San Diego, 
CA, USA). Flow during the analysis was controlled at 
approximately 500 spermatozoa/sec, and 10,000 cells 
were analyzed for each sample. The percentage of FITC- 
positive cells (FL1 channel) was calculated as the per- 
centage of cells with a fluorescence intensity exceeding 
the threshold obtained with the negative control. 

Plasmid construction 

To evaluate the potential effects of PMS2 rsl059060 
(Ser775Asn) polymorphisms on the interaction between 
MLH1 and PMS2, fluorescence resonance energy trans- 
fer (FRET) technology and immunoprecipitation were 
performed. The cDNA encoding MLH1 or PMS2 was 
generated by PCR from a human testis cDNA library. 

For the FRET assay, the primers used for amplifying 
PMS2 (amino acids 655-856) were 5'-CGTT AAGCTT G- 
GAGAAAATCAAGCAGCCGAAG-SVS'-ATACG; 
GATCC CAGGTTGGCGATGTGTCTCAT -3', 
including HindHI and BamHI restriction sites (under- 
lined sequences). Point mutations for PMS2 were per- 
formed using QuikChange Site-Directed Mutagenesis 
Kit (Stratagene, La. Jolla, CA, USA). The amplified frag- 
ment of PMS2 and its genetic variants were cloned into 
the pEYFP-Cl vector (Clonetech, Palo Alto, CA, USA). 
Similarly, the cDNA sequence encoding MLH1 (amino 
acids 506-756) was amplified by PCR using the following 
primers: 5'-CGTT GAATTC GTGTTTTGAGTCTCCAG- 
GAAGAAA-375'-ATAC GGATCC ACACCTCTCAAA- 
GACTTTGTAT-3', which contain EcoRI and BamHI 
restriction sites (underlined sequences). This amplified 
fragment was ligated into pECYP-Cl vector (Clonetech, 
USA). For immunoprecipitation, the cloning of the full- 
length PMS2 and MLH1 cDNA constructs into 
pcDNA3.1 (Invitrogen, Carlsbad, CA, USA), between 
Nhel and BamHI, has already been described [25]. The 
integrity of the inserts was confirmed by sequence 
analysis. 

Cell culture and transfection 

MutLa-deficient HEK293T cells were cultured in 
DMEM: F12 (1:1) (Gibco, Carlsbad, CA, USA), supple- 
mented with 10% foetal bovine serum and 0.1% strepto- 
mycin/penicillin (Gibco, USA) in a humidified 
atmosphere with 5% CO2 at 37°C. Cells were seeded 
onto 30 mm dishes with poly-L-lysine-coated glass cov- 
erslips and co-transfected with YFP recombinant plas- 
mid (YFP-PMS2 or variants of YFP-PMS2) and CFP 
recombinant plasmid (CFP-MLH1) using Lipofectamine 
2000 (Invitrogen) until the cells were at 50 to 60% con- 
fluence, according to the manufacturer's protocols. The 
transfection efficiency was compared by Western blot- 
ting at 72 hours after transfection using anti-PMS2 



(A16-4) (1:100; BD Biosciences), anti-MLHl (G168-728) 
(1:100; BD Biosciences), and anti-p* actin (1:5000; Santa 
Cruz Biotechnology, CA, USA) antibodies. 

Image analysis and calculation of fluorescence resonance 
energy transfer ratios 

We used a Zeiss LSM710 confocal microscope (Carl 
Zeiss, Jena, Germany) operating with a 40 mW argon 
laser. Filter-cube specifications for the fluorescent chan- 
nels were as follows for excitation and emission, respec- 
tively: enhanced cyan fluorescent protein (ECFP), 430 ± 
25 and 470 ± 30 nm; enhanced yellow fluorescent pro- 
tein (EYFP), 500 ± 20 and 535 ± 30 nm; and fluores- 
cence resonance energy transfer (FRET), 430 ± 25 and 
535 ± 30 nm. 

Image analysis involved three basic operations: sub- 
traction of background autofluorescence and blurred 
light, quantification of fluorescence intensity, and calcu- 
lation of a corrected FRET (FRETc) by the following 
equation: 

FRETc = (Ida - a Iaa - d Idd)/Iaa> where I DA is the 
fluorescence intensity from the FRET filter set and I DD 
and I AA are the fluorescent intensities from ECFP (the 
donor) and EYFP (the acceptor), respectively. The cross- 
talk coefficients a and d were considered constant. The 
corrected FRET ratio was defined as FRETc/I DD . 

Co-lmmunoprecipitation and Western blotting 

Proteins were extracted from co-transfected HEK293T 
cells by the M-PER® Mammalian Protein Extraction 
Reagent (Pierce Bio, Thermo, Rockford, IL, USA) 
according to the manufacturer's instruction. Approxi- 
mately 200 ug total cell protein was transferred to a 1.5 
ml microcentrifuge tube, and 20 ul of Protein A/G 
PLUS-Agarose (Santa Cruz Biotechnology, CA, USA) 
was added to the supernatant and the mixture was incu- 
bated at 4°C on a rocker platform for one hour. After 
this incubation, 2 ug anti-MLHl N-20 (Santa Cruz Bio- 
technology, CA, USA) was added and incubated with 
shaking at 4°C overnight. The immunoprecipitates were 
collected by centrifugation at 1,000 x g for 5 minutes at 
4°C, washed 4 times with lysis buffer and then the preci- 
pitates were collected for the Western blotting detection 
with the anti-PMS2 (A16-4) (1:100; BD Biosciences) 
antibody. Proteins were then detected with a Phototope- 
HRP Western Blot Detection kit (Cell Signalling Tech- 
nology, Inc., Beverly, MA, USA). 

Statistical analyses 

Differences in select demographic variables, as well as 
smoking and alcohol status, between the cases and the 
controls, were evaluated using the j 2 test. The Student's 
t test was used to evaluate continuous variables, includ- 
ing age and pack-years of cigarette smoking. The 
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Hardy- Weinberg equilibrium was tested using a good- 
ness-of-fit % test. We used unconditional multivariate 
logistic regression analysis to examine associations 
between genetic polymorphisms and male infertility risk 
by estimating ORs and 95% confidence intervals (95% 
CI). To reduce the potential for spurious findings due to 
multiple testing, we applied the False Discovery Rate 
(FDR) method to the P-values for the differences of gen- 
otype distributions among cases and controls. False Dis- 
covery Rate (FDR) is a new approach to the multiple 
comparisons problem. Instead of controlling the chance 
of any false positives (as Bonferroni methods do), FDR 
controls the expected proportion of false positives 
among suprathreshold voxels [26]. 

Sperm DNA fragmentation was normalized by natural 
logarithm (In) transformation. Linear regression models 
were used to estimate the association with ln-trans- 
formed sperm fragmentation values for each SNP inde- 
pendently. Models were adjusted for age, smoking 
status, drinking status and abstinence time. All P-values 
presented are two-sided and all analyses were carried by 
the Statistical Analysis Software, version 9.1.3 (SAS 
Institute, Cary, NC, USA). 

Results 

Subject characteristics 

The final study population consisted of 1,772 Han Chinese 
individuals, composed of 480 fertile controls, 268 infertility 
patients with non-obstructive azoospermia, 256 infertility 
patients with oligozoospermia (sperm counts < 20 x 10 6 / 
ml) and 768 infertility patients with normal sperm count 



(sperm counts > 20 x 10 6 /ml). The frequency distributions 
of selected characteristics of the case patients and control 
subjects are presented in Table 1. No significant differ- 
ences were observed between cases and controls with 
regard to drinking status or age. However, there was a sig- 
nificantly higher percentage of smokers among cases than 
controls (P < 0.001). Among smokers, cases also reported 
greater cigarette consumption than controls, as assessed 
by the mean number of pack-years (P < 0.05). As expected, 
semen parameters, such as sperm concentration and 
sperm motility, were significantly higher in fertile controls 
than infertile cases (P < 0.001). 

Allelic frequencies and genotype distributions of MMR 
polymorphisms 

The position and minor allele frequency among Chinese of 
the 21 SNPs in the HapMap database are presented in 
Additional file 2 (Table S2). All SNPs were in Hardy- 
Weinberg equilibrium among the controls, except for 
rs3117572 (P = 0.024). Inspection of the cluster plots indi- 
cated good discrimination between genotypes, suggesting 
that these deviations from HWE are likely to be chance 
observations. The genotype distributions among cases and 
controls are presented in Table 2. Overall, the genotype 
frequencies of three SNPs were significantly different 
between the patients with azoospermia or oligozoospermia 
and the controls (P = 0.032 for rs4647269, P = 0.003 for 
rsl059060 and P = 0.002 for rs2075789). Moreover, the 
genotype frequencies of rs 1059060 were also significantly 
different between the patients with normal sperm count 
and the controls (P = 2.0 x 10' 4 ). 



Table 1 Distribution of selected characteristics between cases and fertile controls 


Variables 


Controls 
(n = 480) 


Case 1 a 
(n = 524) 




Case 2 b 
(n = 768) 






N (%) 


N (%) 




N (%) 




Age (mean ± SEM) 


28.1 ± 0.16 


28.3 ±0.16 


0.374 


28.3 ±0.14 


0.355 


Smoking stauts 












Never 


278 (57.9) 


238 (45.4) 


< 0.001 


363 (47.3) 


< 0.001 


Ever 


202 (42.1) 


286 (54.6) 




405 (52.7) 




Pack-years (mean ± SEM) C 


4.3 ± 0.21 


5.2 ± 0.20 


0.002 


4.9 ± 0.14 


0.014 


Drinking status 












Never 


425 (88.5) 


447 (85.3) 


0.129 


667 (86.8) 


0.377 


Ever 


55 (11.5) 


77 (14.7) 




101 (13.2) 




Semen parameters (mean ± SEM) 












Concentration (x 10 6 /ml) 


102.6 ± 3.07 


5.12 ± 0.38 


< 0.001 


73.6 ± 2.12 


< 0.001 


Motility (%) 


65.3 ± 0.58 


3.26 ± 0.27 


< 0.001 


37.9 ± 0.55 


< 0.001 


Volume (ml) 


2.80 ± 0.07 


2.37 ± 0.07 


< 0.001 


2.78 ± 0.05 


0.818 


Sperm DNA fragmentation (%) 


n. d. 


n. d. 




19.5 ± 0.82 





a Case 1: idiopathic infertile men with azoospermia or oligozoospermia. 
b Case 2: idiopathic infertile men with normal sperm count. 
c Among ever smokers. 

d P-values were derived from the x. 2 test for categorical variables (smoking and drinking status) and t test for continuous variables (age and pack-years), 
n. d.: not detected. 
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Table 2 Distribution of the genotype in selected SNPs of 
MMR genes 





tSNP 


f~nntrnl^ 
MAF 


Case1 a 
MAF 




Case2 b 
MAF 




h At 141 

IVlLn 1 
















,. i -70007-7 
rs 1 /yyy// 


0.023 


0.028 


0.628 


0.025 


0.828 




rs4o4/ zoy 


0.047 


0.075 


0.032 


0.067 


0.1 75 




rs 1 54U354 


0.304 


0.343 


0.326 


0.33 1 


0.428 


rlvlbz 
















rsio 1 jjoj 


0.332 


0.332 


0.923 


0.342 


0.81 6 




rszzooooU 


n 07y1 

U.U/4 


U.UoU 


0 111 
U./ /z 


n oqo 


u.4y i 




rs 1 1 /by^ou 


0.41 0 


0.41 9 


0.81 6 


0.379 


0.736 




rr- 1 OCOO^O 

rs i loyuou 


0.091 


0.1 46 


0.003 


0.1 70 


2.0 x 
10" 4 




rs2228006 


0.063 


0.054 


0.520 


0.056 


0.582 


MLH3 
















rsl 371 2 


0.184 


0.185 


0.913 


0.185 


0.962 




rs71 56586 


0.222 


0.217 


0.832 


0.406 


0.527 




rsl 75049 


0.185 


0.173 


0.646 


0.177 


0.727 


MSH4 
















rs 102 1462 


0.309 


0.336 


0.592 


0.317 


0.842 


MSH5 
















rs3749953 


0.121 


0.154 


0.278 


0.148 


0.653 




rsl 150793 


0.142 


0.141 


0.960 


0.138 


0.827 




rs707939 


0.359 


0.337 


0.582 


0.357 


0.886 




rs707938 


0.304 


0.331 


0.557 


0.336 


0.231 




rs31 15672 


0.375 


0.402 


0.681 


0.377 


0.929 




rs31 17572 


0.223 


0.206 


0.724 


0.231 


0.782 




rs2299850 


0.063 


0.056 


0.652 


0.068 


0.720 




rs9461718 


0.144 


0.161 


0.472 


0.157 


0.567 




rs2075789 


0.081 


0.139 


0.002 


0.097 


0.485 



Abbreviations: MAF, minor allele frequency. 

a Case 1: idiopathic infertile men with azoospermia or oligozoospermia. 
b Case 2: idiopathic infertile men with normal sperm count. 
Data in bold highlights the statistic significant results. 
c False Discovery Rate (FDR) corrected P-value. 
Data in boldface represent P < 0.05. 

Logistic regression analyses showed that in the domi- 
nant-effect model, significantly increased risks of azoos- 
permia or oligozoospermia were associated with 
rs4647269 CT/TT (adjusted OR = 1.63, 95% CI: 1.10 to 
2.41), rsl059060 GA/AA (adjusted OR = 1.60, 95% CI: 
1.17 to 2.18) and rs2075789GA/AA (adjusted OR = 
1.83, 95% CI: 1.32 to 2.55), as compared to wild-type 
homozygous carriers (Table 3). Meanwhile, a signifi- 
cantly increased risk of male infertility with normal 
sperm count was associated with the rsl059060 GA/AA 
genotypes (adjusted OR = 1.83, 95% CI: 1.37 to 2.43), as 
compared to wild-type homozygotes. 

Association between MMR polymorphisms and sperm 
DNA fragmentation 

Considering the importance of the MMR pathway in 
maintenance of DNA integrity, we further evaluated the 



effects of these three SNPs on sperm DNA fragmenta- 
tion. In the present study, semen samples were pre-trea- 
ted with cryopreservation prior to TUNEL analyses. 
However, it has been demonstrated that the process of 
cryopreservation can lead to an increase in oxidative 
stress and percentage DNA fragmentation [27]. To 
determine whether the results of the TUNEL analyses 
were profoundly influenced by cryopreservation in our 
study, 10 semen samples were pre-treated with or with- 
out cryopreservation prior to TUNEL analyses. As 
shown in Additional file 3 (Table S3), modest but sig- 
nificant elevated levels of sperm DNA fragmentation 
were induced by cryopreservation (P = 0.001). However, 
all the semen samples undergo the same cryopreserva- 
tion process, thus we believe that the effect of cryopre- 
servation, if any, is unlikely to be substantial. The non- 
normal distribution of sperm DNA damage levels and 
sperm concentration were natural log (ln) transformed 
for further association studies (skewness-kurtosis tests P 
> 0.05). After adjustment for age, smoking, alcohol use 
and length of abstinence, we found that subjects who 
carried the rs4647269 CT/TT genotypes displayed 
markedly higher levels of sperm DNA fragmentation 
compared with the CC homozygotes (mean ± S.D., 
11.82% ± 2.66% vs. 26.58% ± 1.97%; P < 0.001) (Figure 
1A). Moreover, a gradual increase in sperm DNA frag- 
mentation was found among the three PMS2 rsl059060 
subgroups (mean ± S.D., 12.30% ± 2.72%, 17.99% ± 
2.27%, and 24.78% ± 1.70% for GG, GA and AA, 
respectively; P tre nd < 0.001) (Figure IB). However, no 
significant difference was observed for the MSH5 
rs2075789 (Figure 1C). 

Effects of the PMS2 Ser775Asn polymorphism on MLH1 
and PMS2 interaction 

The PMS2 Ser775Asn polymorphism (rsl059060) was 
potentially located within the MLH1-PMS2 interacting 
domain. Therefore, we examined whether PMS2 Ser775- 
Asn polymorphisms influence binding between MLH1 
and PMS2. HEK293T cells were transiently co-trans- 
fected with plasmids encoding the MLH1 (amino acids 
506-756) and wild-type or genetic variants of PMS2 
(amino acids 675-850). The schematic diagram of the 
FRET assay is summarized in Additional file 4 (Figure 
SI). By confocal fluorescence detection, we found that 
there was a weak interaction between MLH1 and 
PMS2-775Asn proteins, for little FRETc was detected in 
cells co-expressing CFP-MLH1 and YFP-PMS2-775Asn 
plasmid (Figure 2B). Cells co-expressing CFP-MLH1 + 
YFP-PMS2-775Ser had a four-fold increase in FRETc 
values (0.031 ± 0.013, n = 12; 0.008 ± 0.005, n = 12; P < 
0.001) (Figure 2A). This result suggested that the PMS2 
Ser775Asn polymorphism could significantly influence 
the interaction between MLH1 and PMS2. 
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Table 3 Genotype frequencies of the four SNPs in MMR genes in patients and controls and their associations with 
male infertility risk 



SNP ID 


Genotype 


Controls (n = 480) 


Case 1 a 


(n = 524) 


Case 2 b 


(n = 768) 


N (%) 


N (%) 


OR (95% Cl) c 


N (%) 


OR (95% Cl) c 


MLH1 














rs4647269 


CC 


431 (905) 


444 (85.5) 


Reference 


665 (86.6) 


Reference 




CT 


45 (94) 


72 (13.9) 


1 .56 (1 .05 to 2.32) 


94 (12.2) 


1 .34 (0.92 to 1 .96) 




TT 


0 (0.0) 


3 (0.6) 


NA 


4 (0.5) 


NA 




CT/TT 


45 (9.4) 


75 (14.4) 


1.63 (1.10 to 2.41) 


98 (12.7) 


1.39 (0.93 to 2.01) 


Ptrend 








0.009 




0.044 


PMS2 














rs 1059060 


GG 


393 (82.6) 


387 (74.6) 


Reference 


534 (69.5) 


Reference 




GA 


79 (16.6) 


112 (21.6) 


1.43 (1.03 to 1.96) 


197 (25.6) 


1.82 (1.36 to 2.44) 




AA 


4 (0.8) 


20 (3.8) 
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P trend 
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0.0002 




0.224 



a Case 1: idiopathic infertile men with azoospermia or oligozoospermia. 
b Case 2: idiopathic infertile men with normal sperm count. 
Adjustment for age, smoking status and alcohol use. 



We also used a co-immunoprecipitation assay to detect 
the effects of PMS2 variants on the MLH1 and PMS2 
interaction. Full-length MLH1 and PMS2 775Ser or 
PMS2 775Asn were constructed and co-transfected into 
HEK293T cells. Western blot analysis of whole cell 



lysates showed satisfied transfection efficiency (Figure 
3A). The co-immunoprecipitated result also suggested 
that binding between MLH1 and PMS2-775Ser was 
more robust compared with binding between MLH1 
and the PMS2-775Asn variant (Figure 3B, lane 4 vs. 3) 
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Figure 1 Box-and-whisker plots of sperm DNA fragmentation for different genotypes The boxes represent the 25 th and 75 th percentiles; 
whiskers are lines extending from each end of the box covering the extent of the data on 1.5 x inter-quartile range. The median value is 
denoted as the line that bisects the boxes. Circles and asterisks represent the outlier values. Significant differences were measured by multiple 
inear regression. 
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Figure 2 FRET imaging of MLH1 and PMS2 interaction in live HEK293T cells Images of CFP-tagged (green) and YFP-tagged (red) 
constructs when transiently expressed in HEK293T cells. Co-localization of co-expressed constructs is shown as yellow in overlay images. The 
pseudocoloured images represent FRET signals corrected for any bleed-through using the micro-FRET method (FRETc). A: Co-localization 
(overlay) and direct interactions (FRETc) between MLH1-CFP + PMS2 (wt)-YFP were detected in the nucleus. B: Cells co-expressing MLH1-CFP + 
PMS2 S775N-YFP showed good co-localization of fluorescent signals but little detectable FRETc signal in the nucleus. 

* y 



which suggested that mutation of PMS2 significantly 
attenuated the protein-protein interaction of MLH1 and 
PMS2. 

Discussion 

Accumulating evidence demonstrates that MMR plays a 
critical role in the maintenance of genetic integrity and 
participates in the meiotic recombination process 
[14-16]. Although mutations in MMR genes are 



considered as potential risk factors for various cancers 
[28,29], only limited data exist on the potential role of 
polymorphisms in the MMR genes on male infertility. 
To our knowledge, this study is the first to provide a 
comprehensive evaluation of the relationship between 
polymorphisms in MMR and susceptibility to male 
infertility in a relatively large sample size. On the basis 
of analysis of 480 controls and 524 infertility patients 
with azoospermia or oligozoospermia, we observed that 




Figure 3 Interaction studies between hMLH1 and hPMS2 variants. A: Western blot of total protein extracts (50 ug each) from HEK293T cells 

transfected with pcDNA3.1-MLH1 and either wild-type pcDNA3.1-PMS2 (775Ser) or pcDNA3.1-PMS2 (775Asn) variants, p-actin was used as 

controls. B: The lysates of cells co-expressing the two plasmid were immunoprecipitated with anti-MLH1 N-20 antibody, and then detected with 

anti-PMS2 (A16-4) antibody. Western blot signals were quantified employing Quantity-One software. 
\ . J 
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one intronic SNP in MLH1 (rs4647269) and two non- 
synonymous SNPs in PMS2 (rsl059060, Ser775Asn) and 
MSH5 (rs2075789, Pro29Ser) were associated with 
increased susceptibility to poor sperm production. 

As an important pathway in the DNA damage repair 
network, MMR also plays a critical role in the mainte- 
nance of genetic integrity. Thus, it would be expected 
that these three significant SNPs that affect sperm DNA 
integrity could also modify male infertility risk. Based on 
a case-control study consisting of 480 controls and 768 
patients with normal sperm count, we found that PMS2 
rsl059060 was significantly associated with male inferti- 
lity with normal sperm count. Further analysis based on 
450 infertile men revealed significant associations of 
MLH1 rs4647269 and PMS2 rsl059060 with sperm 
DNA fragmentation. However, we did not detect any 
association between MSH5 Pro29Ser polymorphisms 
and sperm DNA damage. This result is explained by the 
fact that MSH5 is a meiosis-specific protein crucial for 
reciprocal recombination, and it has no apparent mis- 
match repair activity because it is missing the appropri- 
ate amino acid residues [30]. 

MLH1 and PMS2 form the MutLot heterodimer that 
leads to the repair of mismatched DNA through activa- 
tion of exonuclease-mediated degradation of DNA [31]. 
Guerrette et al. localized the MLH1-PMS2 interaction 
region to amino acids 506-675 of MLH1 and amino 
acids 675-850 of PMS2 [32]. It is conceivable that the 
PMS2 Ser775Asn polymorphism could directly impact 
the integrity of the interaction between MLH1 and 
PMS2. In the present study, we provided evidence, for 
the first time, that the PMS2 Ser775Asn variant attenu- 
ates the interaction of MLH1 and PMS2, as illustrated 
by FRET and co-immunoprecipitation assays. 

The MSHS rs2075789 polymorphism in the coding 
region of the human MSHS gene leads to a proline to 
serine alteration and is located within the MSH4-MSH5 
interacting domain. To address the effect of the Pro29- 
Ser alteration on the interaction between MSH4 and 
MSH5, a quantitative yeast two-hybrid analysis has been 
performed [33]. This alteration causes a moderate but 
significant reduction in the interactions between both 
proteins, which could affect the formation of the 
MSH4-MSH5 heterocomplex. These findings strongly 
support our molecular epidemiological observation that 
the MSH5 Pro29Ser polymorphism is associated with a 
significantly increased risk of azoospermia or oligozoos- 
permia. Supporting evidence also comes from associa- 
tion studies by other investigators. In a recent study of a 
Chinese population with a small sample size, Xu et al. 
observed a 2.89-fold increased risk of azoospermia or 
oligozoospermia among the MSHS Pro29Ser allele car- 
riers [34]. In addition, a case-control study including 41 
women with premature ovarian failure and 39 controls 



suggested that there is a correlation between the MSH5 
Pro29Ser polymorphism and premature ovarian failure 
in women [35]. 

Another SNP associated with risk in our study 
(rs4647269) is intronic. However, SNP rs4647269 tags 
SNP rs9852810 (r 2 = 1, D' = 1), which was associated 
with prostate cancer risk and prostate cancer recurrence 
[36]. Because both of these two SNPs are located in the 
intron of the MLH1 gene, it is uncertain which one of 
these two variants causes increases in male infertility 
risk. To identify additional SNPs that could be asso- 
ciated with male infertility risk that may be in high link- 
age disequilibrium (LD) with these two sites, we 
screened all of the common variants (with MAF > 0.05) 
within an approximately 20 kb-long region surrounding 
these two sites (approximately 10 kb upstream and 
approximately 10 kb downstream of these loci) based on 
the CHB HapMap data resource. We found that 
rs4647269 is in complete LD with SNP rsl046512, 
which is located approximately 2.5 kb upstream of start 
codon of MLH1. Therefore, it is highly likely that the 
rsl046512 SNP near the 5' region of the MLH1 gene 
may be the causal variant. 

Another interesting finding was that smoking was 
associated with increased risk of male infertility. 
Although the effects of tobacco cigarette smoke on male 
reproduction are somewhat inconclusive, a number of 
studies have shown higher incidences of abnormal 
sperm morphology [37,38] and decreased sperm motility 
concentration in men who smoke [39,40]. A meta-analy- 
sis [41], including 27 studies, indicated that cigarette 
smoking is associated with a 13% reduction in sperm 
concentration, a 10% reduction of sperm motility, and a 
3% reduction of morphologically normal sperm. Further- 
more, fluctuation in reproductive hormone levels have 
been documented in male smokers [42,43]. However, 
the mechanism(s) of these changes, if any, remains 
unclear. 

Of note, like all case-control studies, selection bias 
may exist and might influence interpretation of the 
results. However, we believe that potential confounding 
bias might have been minimized by matching the con- 
trols to the cases on age and by further adjustment for 
the confounding factors in statistical analyses. In addi- 
tion, the fact that genotype frequencies of all SNPs in 
our controls fit Hardy- Weinberg equilibrium and were 
similar to those obtained from the HapMap Project 
further supports the randomness of our control selec- 
tion. We believe that the selection bias, if any, is unli- 
kely to be substantial. 

Conclusions 

The present study extends the previous understanding 
of the MMR polymorphisms and their effects on the 
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risk of idiopathic azoospermia or oligozoospermia by 
further evaluating the contribution of these polymorph- 
isms in relation to sperm DNA fragmentation. These 
novel findings might be helpful in improving the under- 
standing of the role of genetic variation in susceptibility 
to reduced sperm DNA integrity and in providing diag- 
nostic implications for assisted reproduction success 
rates. Although, these three SNPs (rs4647269, rsl059060 
and rs2075789) associated with risk in our study are sig- 
nificantly higher for some variants in the patient group, 
the actual rates are quite low and would potentially 
account for a low percentage of infertility. However, it is 
important to know that genetic variants associated with 
common complex diseases like male infertility are only 
"one piece of the puzzle" making up an individual's 
overall risk for disease. It is highly likely that the genetic 
risk for developing male infertility is influenced by the 
additive effects of many different genetic variants and 
other risk factors. So, further research is required to 
define their interactions with other susceptibility alleles 
and environmental factors can lead to a substantial 
increase in male infertility risk, especially when exposed 
to certain dietary and lifestyle habits. 

Additional material 
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MMR: mismatch repair; SNP: single-nucleotide polymorphism; TUNEL: Tdt- 
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